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ABSTRACT 

We offer a new method for determining the wind source term for energy and momentum fluxes 
transfer from the atmosphere to the wind-driven sea. This new source-term formulation is 
based on extensive analysis of experimental data collected at different sites around the world. 
It is shown that this new wind source term to be consistent both with numerical solution of 
exact equation for resonant four-wave interactions and available experimental data. 
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1. Introduction 



Presently, most oceanographers believe that a sat isfactory d e script i on of the the ey o 



of wave action in a wind -driven sea is given by (see 



and 



Badulin et al. 



(120071 1): 



Zakharovl fcoOSf ) 



Badulin et al. 



ution 



fcoosh 



dNi, dNk duk 9iVk 



dt 



dt 



9k dr 



(1) 



where = A^(k, r, t) is the wave action density, k = {kx,ky) is the Fourier space wave- 
number vector, r = {x, y) is the real-space coordinate vector, ojk = y/gk (where k = |k|) is the 
dispersion law for gravity waves, Sni is the source term for resonant four-wave interactions, 
and Ss = Sin + Sdiss is the source term responsible for energy transfer from the wind and 
dissipation due to wave breaking and other forms of dissipation within the wave field. 

The Sni term is a nonlinear integral operator which, in a global sense, can be written as 
a conservative kinetic equation 



dN^ 
dt 



5, 



nl 



(2) 



This form conserves total wave action and formally conserves energ y and momentum; how- 



ever, c onservation is not realized over a finite range of frequencies (see lPushkarev and Zakharov 



(j2000[ )). The associated nonlinear fluxes have been shown to produce power Kolmogorov- 
like tails in a finite time, transferring energy and momentum into the area of small wave 
numbers. The first spa t ially homogeneous isotropic solution of that type was found by 



Zakharov and Filonenkd (119661 ) 



pl/3 

r' 



(3) 



where spectral energy density e{u, 6) {6 is the wave-number vector angle) is connected with 
wave action density iV(k, t) by the relation 



e{uj, 9)dud9 = UkN^dk 



(4) 



Pq is the energy flu x to high wave numbers and Ci is the Kolmogorov constant Ci 



47r ■ 0.219 = 2.75 (see 



Zakharovl tom ). 



In contrast to the strong theoretical basis for Sni, our knowledge of the wind input and 
wave dissipation source terms Sg is quite poor. The creation of a reliable, well justifled 
theory for wind input SiU has been hindered by the presence of strong turbulence in the 
air boundary layer over the sea surface. Even the one of the most crucial elements for this 
theory, the vertical proflle of the mean horizontal wind velocity, is still subject to debate in 
the region closest to the ocean surface where wave motions interact strongly with atmospheric 
motions. 

Direct measurements of Sm are scarce. As a result, we currently have many different 
heuristic models of and the scatter between different models is large in terms of both the 
predicted tot al fluxes and the distribution o f fluxes into different frequencies and directions. 



For instance, 



Donelan and Pierson-ir 



imately flve times higher then the 



19871) form for Sin pre dicts magnitudes of Sm approx- 



Hsiao and Shemdin 



different forms for of 5'^ is presented by 



Badulin et al 



983f) form for Sin- A comparison of 



fl2005h . 



Unfortunately, but our understanding and quantiflcation of the dissipation source term 
Sdiss is not any better. The theory is not well developed and experimental data are far from 
complete; consequently, the forms for Sdiss used in operational models are heuristic and not 
well justified. 



The situation in wave models is aggravated by the fact that the heuristic forms for 
wind input and wave dissipation are of the same magnitude as the nonhnear source term in 
V all parts of the wave spectrum, i.e. Sg ~ Sni- This viewpoint was postulated by 
and s in ce this time has been commonly accepted. Some authors, for example. 



essentia] 


V all 


Phillips 


(1985) 


Donelan et al. 



(119851) 



Donelan] ( 119991 )) even assume that Sni is negligibly small with respect 
to Ss- According to this opinion, the wind-generation process at full development is arrested 
by wave-breaking; although it is difficult to find strong theoretical or empirical support for 
these hypotheses. 

Although the magnitude of Sni can be small in regions of a spectrum exhibiting constant 
fiuxes, in many situations the nonlinear term Sni surpasses the Sg source term b y at least 



the or d er of magnitude in many regi ons of the spectrum. It has been argued by IZakharov 



torn 



Zakharov and Badulid (120111 ) that this makes the role of Sni absolutely critical to 
the evolution of wind wave spectra. In the next three sections, we shall utilize experimental 
evidence, theoretical considerations and numerical simulations to formulate and test a new 
wind input source term that functions in concert with Sni to produce wave growth and wave 
spectra in agreement with the findings of Badulin et al. (2005, 2007). 



2. Experimental evidence 



We start here by examining empirical evidence from around the world whi ch has been 



utilize d to quantify energy levels within the equilibrium range in spe ctra by 



(120041 ). For convenience, we shall also use the same notation used by 



Resio et a" 



Resio et al 



(I2004h 



in their study, for the directionally integrated spectral energy densities in frequency and 



wavenumber space, 



2iia4Vg 



(5) 
(6) 



where f = is the constant, V is some characteristic velocity and (3 = ^a^Vg-^^'^. 

These notations are based on relation of spectral densities E{f) and F{k) in frequency 
/ = ^ and wave-number k basises: 

F{k) = ^E{f) (7) 

where Cg = ^ is the group velocity. 

The notations in Eqs.(|5])-([6]) are connected with the spectral energy density e{uj,6) (see 
Eqs.dSJ-dl])) through 



E{f) = 27ry" e{uj,e)de 



(8) 



The iResio et al.l ( 120041 ) analysis showed that experimental energy spectra F{k), estimated 
through averaging < k'/''F{k) >, can be approximated by linear regression line as a function 
of {u\cpY/^g~^/''. FigH] shows that the regression line 



1 



(3 - -"4 [("A'-P 



Mo] g 



-1/2 



(9) 



indeed, seems to be a reasonable approximation of these observations. Here 04 = 0.00553, 
Mo = 1.93 m/ sec, Cp is the spectral peak phase speed and ux is the wind speed at the elevation 
equal to a fixed fraction A = 0.065 of the spectr al peak wavelength 2iT/kp, where kp is the 



spectral peak wave number. 



Resio et al. 



(j2004| l assumed that the near surface boundary 



layer can be treated as neutral and thus follow a conventional logarithmic profile 



having Von Karman coefficient k = 0.41, where z = X- 2'K/kp is the elevation equal to a fixed 
fraction A = 0.065 of the spectral pea k waveleng t h 27r/ kp, where kp is the spectral peak wave 



number, and zq = acul/g subject to 



Charnockl (119551 ) surface roughness with ac = 0.015. 



3. Theoretical considerations 



S elf- similar solut ions of conservative kineti c equation ([2]) were studied in lZakharovl (j2005[ ) 



Baduhn et al. 



(120051) and 



Badulin et al. 



( 120071 ). In this chapter we study self-similar solutions 



of the forced kinetic equation 



de{oj, 



dt 



Sni + j{uj,9)e{u,9) 



(11) 



where e(w, 6) = ^N(k, 6) is the energy spectrum. For our purposes, it is sufficient to simply 



use the dimensional estimate for S, 



nil 



Eq. (fTT]) has a self-similar solution if 



(12) 



(13) 



where s is a constant. Looking for self-similar solution in the form 



e{u,t) = tP+''F{ujt'') 



(14) 



we find 



P 



1 



s + 1 

9q-l 8-s 



2 2{s + l] 
6 



(15) 
(16) 



The function F{C,) has the maximum at ^ ~ ^p, thus the frequency of the spectral peak is 

oop ^ ^pr" (17) 

The phase velocity at the spectral peak is 

= JL = ^i, = f (18) 

According to experimental data, the main energy input into the spectrum occurs in the 
vicinity of the spectral peak, i.e. at w ~ Up. For u » Up, the spectrum is described by 
Zakharov-Filonenko tail 

e{oo) ~ P'/^io~^ (19) 



Here 



P = I / j{uj,9)e{Lo,9)de (20) 



'0 JO 

This integral converge if s < 2. For large u 



e(w,^)-^-^^ (21) 



More accurately 

e(a;,t) ^ ^«^-''c>(^) (22) 
2-s 

V = (23) 

Now supposing s = 4/3 and 7 ~ cu^/'^, we get 77 = 1/3 , which is exactly experimental 
regression line prediction. Because it is known from regression line on Fig{T]that ^ = 1/3, 
we immediately get s = 4/3 and the wind input term 

S^,nd ^ cu'/'e (24) 



One should note that the dependence shown in Eq. flM|) has already been predicted by 
Resio and Perrid (119891 ) from dimensional consideration. 



4. Numerical simulations 

To check the self-similar conjecture implicit in our formulation shown in Eq.f l24p we 
performed a series of numerical simulations of Eq.([T]) in the time domain = 0, or 
spatially homogeneous) and fetch limited = 0, or spatially inhomogeneous) situations. 
The same input term Eq. (1241) has been used in both cases in the form 

4/3 



wind 



f{0) 



0.2 



Pair 


( UJ \ 






Pwater 




COS 9 


for 



-n/2 <9 <n/2 



otherwise 
1.3- 10" 



9 Pair -I Q -I 



^ Pwater 

where u is the wind speed, pair and plater are the air and water density correspondingly. 
Wind speed u is taken here as the spe ed at a refer e nce le vel of 10 meters (uio)- To make 



com parison w i th ex perimental results of 



see 



Resio et al. 



( 120041 ). we used relation ~ "^^10/28 



Gohtsvnl mm ) in Eq.([l0D. 
Both situations also need knowledge of the dissipation term Sdiss, which is taken into 
account here in the same way, as was proposed by Resio and Long (2007), where white- 
capping dissipation term Sdiss was introduced implicitly through an energy spectral tail 
stretching in frequency range from = 1.1 to fmax = 2.0. To date, this a pproach has 



been shown by both experimental observations and numerical experiments (see 



Resio et al. 



(|2004| )) to provide an effective sink for energy in direct cascades into high frequencies. 



a. Time domain simulations 



Time domain numerical simulations have been started from an initial uniform noise 
energy distribution in Fourier space. The typical picture of the directional energy spectrum 
for the wind speed u = 10m/ sec is presented on Figj2j In this figure, one can distinguish 
three separate areas of the spectrum area of spectral peak, intermediate portion of the 
spectral tail proportional to and high-frequency portion of the spectrum, proportional 
to 

The compensated spectrum F{k) ■ k^^'^ from these simulations is presented on Figj3l One 



can see plateu-like region responsible for k behavior, equiva| 



exact solution of Eq.([T]) was found by iZakharov and Filonenko 



e nt to /-^ tail in FigJl This 



fll966h 



Now we are tu rning to direct com parison of the numerical simulation with the experi- 



mental analysis by 



Resio et al. 



(j2004j ) presented on FigJH 
FiglH presents the plot of the function (3 = F{k) ■ k^/'^ in terms of {u\CpY^'^ / g^^"^ for four 
different runs, corre sponding to wind speeds u = 2.5,5.0,10.0,20.0 m/sec, along with the 



regression line from 



Resio et al. 



(120041 ). The overall correspondence with the regression line 
for values of n = 2.5,5.0, 10.0 m/sec is quite good. The results corresponding to n = 20.0 
m/sec are bit off of the regression line, but exhibit the same slope. 

Another important theoretical relationship can be derived from joint consideration of 
Eqs. dSD, dZD and ([22]): 



1000/3 = 3 



V ^1/3 



[u^Cp) 



-.1/2 



(25) 



This dependence, presented on FigJU shows that using this new wind input term produces 
a good correspondence of theory, experiment and numerical results. 



b. Limited fetch numerical simulations 



A limited fetch equivalence to the time-domain situation can be described by reduction 



where x is the coordinate orthogonal to the shore and 6 is the angle between individual 
wavenumber k and the axis x. 

Eq.f l26|l is somewhat difficult for nimerical simulation, since it contains a singularity in 
the form of cos in front of H. We overcame this problem of division by zero through zeroing 
one half of the Fourier space of the system for the waves propagating toward the shore. Since 
it is well- known that the energy in such waves is small with respect to waves propagating 
in the offshore direction, such approximation is quite reasonable for our purposes. 

The same sort of self-similar analysis performed for time domain situation can be also 
repeated for limited fetch one. The result of inserting the self-similar substitution 



of Eq.© 




(26) 



(27) 



into Eq. fl26|) gives the values of indices 



P = 



lOg- 1 
2 



s + 2 



1 




and 



4 



(29) 



^=10' 
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We find that in the fetch-limited case the wind forcing index s is similar to time domain 



case when the wind forcing is given by Eq. (lM|) . Therefore, numerical simulation of Eq.f l26p 
has been performed for the same input functions as in the time doman case with the same 
initial conditions for the form of low-level energy noise in Fourier space. 

FiglSpresents characteristic spectral energy distribution for our hmited fetch simulations. 
Figj6] presents the plot of function (3 = F{k) ■ k^^"^ in terms of the parameter {u\CpY^^ / g^^"^ 
for two different runs, corresponding to the same wind speed u = 10.0 m/sec: one for time 
domain simulations and another for limited fetch simulations. Both sets of simulations show 



good correspondence with the 



Resio et al. 



( 120041 ) regression line. 



5. Conclusion 

We have inroduced a new form 
considerations by 



'or wo nd input, similar to that derived from dimensional 



Resio and Perrid ( 119891 ). Numerical simulations using this new source term 



show good agreement with observed equilibrium range characteristics around the world and 



produce spectral shapes and energy 



theoretical arguments of 



Zakharovl (l2005f ) 



grow t h rates which are con sisten t with the sel 



Baduhn et al. 



(l2005h and 



Baduhn et al. 



-similar 



(120071 ) 



It is believed that this new form of wind input term can improve the quality of ocean wave 



prediction in operational models forecasts. 
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Fig. 1. Correlation of equilibrium r ange coeffici e nt w ith {uj^CpY^^ / g^^'^ based on data from 
six disparate sources. Adopted from \Resio et al. 
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Fig. 4. Experimental, theoretical and numerical evidence on the single graph for 1000/3 as 
a function of (uj^CpY^^ / g^^"^ . Experimental result: dotted line - experimental regression line 
from FigUl Theoretical result: dashed line - theoretical relation Eg. i[25\). Numerical results: 
crosses correspond to u = 2.5, stars to u = 5.0, rectangles to u = 10.0, triangles to u = 20.0 
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Fig. 5. Energy spectral density e{f,6) as a function of frequency f and angle 9 in polar 
coordinates. 
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Fig. 6. Two simulations for wind input 10m/ sec : time-limited domain (triangles) and lim- 
ited fetch (crosses). Dotted line - correlation of the equilibrium range coefficient (3 with 
{u\cpY^^ / g^^"^ based on data from six disparate sources adopted from Resio at al, 2004- 
Dashed line - theoretical value of equilibrium range coefficient (i. 
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